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I will discuss recent results of the NPLQCD Collaboration regarding the calculation of hadronic 
interactions with lattice QCD. A particular emphasis will be spent on pi-pi scattering and other 
meson interactions. 
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1. Two-Particle Energy Levels in a Finite Euclidean Volume 

In a finite Euclidean volume, one does not calculate scattering matrix elements directly [1], 
but rather one calculates the volume dependent distortion of the two-particle energy levels due to 
the particle interactions which, below inelastic thresholds, can be related to the infinite volume 
scattering phase shift; this is known as Luscher's method [2]. For example, for two degenerate 
particles in a cubic spatial volume of size and periodic boundary conditions, the distortion of the 
free-particle energy eigenvalues is related to the scattering phase shift through through a regulated 
three-dimensional summation, 

r-^ ^ o pCOtdip) 1 ( V^L^\ 

M:2,n=2^p^+m^-lm^ ^ ^ = jS('—^] (1.1) 

where S{z) can be found in Refs. [3], 

Tn^-4^A. (1.2) 
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2. Precision I = lTnz Scattering 

Low energy nn scattering is theoretically the best understood hadron-hadron interaction. At 
leading order in chiral perturbation theory (;t^PT) [4, 5], the effective field theory describing the 
low-energy interactions of QCD, the s-wave scattering lengths are uniquely predicted [4]. The sub- 
leading corrections are known to one- [5] and two-loops [6]. Each of these higher order corrections 
depends upon unknown low-energy-constants (LECs) which must be determined from comparison 
with experiment or lattice calculations. 

Low energy nn interactions are also numerically very tractable as the pions are the lightest 
hadrons of QCD. The exotic I = 2 Tin scattering channel ^ is therefore an ideal process to study with 
lattice QCD before exploring the more challenging pion-nucleon and nucleon-nucleon interactions. 
The NPLQCD Collaboration performed the first lattice calculation of the 7 = 2 nn scattering length 
with 2-1-1 dynamical flavors of light quarks {up, down, and strange) [7]. Using continuum x?T to 
perform the extrapolation to the chai^ged pion mass, they found nij^a^j^ = —0.0426 ± 0.0019. The 
calculation was performed by computing domain-wall valence quark propagators [8] on the Asqtad 
improved [9] coarse MILC lattices [10], the mixed-action lattice scheme of the LHP Collabora- 
tion [11]. Subsequently, a high-precision update of this result has been published [12], finding at 
the chai^ged pion mass 

m;j+ ^ -0.04330 ± 0.00042. (2.1) 

There are several ingredients which were necessary for this precision calculation; high statistics 
with the lattice calculation, systematic treatment of the lattice-artefacts from the discretization and 
finite volume, and chiral symmetry. A series of papers allowed for all of these systematic eiTors 
to be addressed in the updated calculation [12]. First, Ref. [12] included an order of magnitude 
increase in statistics compai^ed to the original calculation [7], see Table 1 for details. In Refs. [13], 
the lattice spacing corrections to the scattering length were determined by utilizing the relevant 

'Due to the disconnected quark diagrams, the scalar / = channel is numerically more expensive and noisier. 
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Table 1: Parameters of NPLQCD lattice calculation of / = 2 nn scattering [12], using coarse MILC ensem- 
bles [10] with lattice spacing b ^ 0.125 fm and volume L ^ 2.5 fm. 



{m^,mK}[yyeV] 

Ncfg X N.ource 



{290,580} {350,595} {490,640} {590,675} 
468 X 16 658 x 20 486 x 24 564 x 8 



mixed-action ;tPT. In Refs. [14], the exponentially suppressed volume corrections to Eq. (1.1) 
were determined.^ The other important ingredient is chiral symmetry, both the chiral symmetry 
of the domain-wall valence propagators used in the computation, as well as the fact that the Tin 
scattering lengths must vanish in the chiral limit. When expressed in terms of the lattice-physical 
pion mass and decay constant, the mixed-action extrapolation formula for the scattering length is 
given at one loop by [13] 

a4 



-{2n) 



(47r/;r) 
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(2.2) 



where l'j^{pL) is the same linear combination of renormalized Gasser-Leutwyler coefficients as 
appears in continuum ;^PT. There is no explicit dependence upon the valence-sea mesons [16] and 
the only explicit sea quark effects are contained in the last term, where for the quark mass tuning 
employed in Refs. [7, 12], Ap^ = b^\\, which is the staggered taste-Identity splitting, which has 
been computed on the coarse MILC lattices [17]. The resulting chiral extrapolation is presented 
in Fig. 1 . One striking feature of this lattice calculation, for the full range of quark masses used, 
the resulting scattering length is consistent Weinberg's tree level prediction [4]. Aside from the 
quantum corrections to the pion mass and decay constant, there is no evidence of higher order 
effects. This feature is not specific to / = 2 nn scattering but has also been found in computations 
of both the I = 2> /I Kn anA I = I KK scattering lengths [18]. 

3. Nuclear and Hyper-Nuclear Interactions 

The main objective of the NPLQCD Collaboration is to make predictions for the structure and 
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Figure 1: Chiral extrapolation of I = 2 mt scattering length using mixed-action XPT. 



^The volume corrections computed in Ref. [14] were determined with continuum ;tPT. The volume corrections for 
partially-quenched/mixed-action theories are also known [15] and were used in Ref. [12]. 
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Figure 2: '5o and ^Si nucleon-nucleon correlation functions taken from Ref. [19]. 

interactions of nuclei using lattice QCD. To that end, the first dynamical calculations of nucleon- 
nucleon and hyperon-nucleon scattering lengths were performed recently [19]. However, lattice 
calculations of (multi)-baryon correlation functions suffer from a exponential degradation of the 
signal-to-noise [20], as can be seen in Fig. 2. This noise degradation is currently the main obstacle 
to computing multi-baryon correlation functions at light pion masses. Solutions to this problem 
will require several orders of magnitude increase in statistics as compared to the calculations in 
Refs. [19], especially as one utilizes lighter quark masses. There has been a proposal to alleviate 
this noise problem with a clever choice of boundary conditions [21], although this proposal has yet 
to be implemented numerically. 

4. Kaon Condensation 

In addition to calculating properties of two pions or two nucleons, the NPLQCD Collaboration 
has utilized these same computing resources to calculate properties of up to twelve pions and 
kaons [22, 23]. In order to analyze these correlation functions, the energy levels of n non-relativistic 
identical bosons interacting in a periodic cubic spatial volume of length L was determined through 
ff{L~'^) [24]. For example, we were able to compute the kaon chemical potential as a function the 
density of K^s at densities relevant for neutron star phenomenology. Additionally, three -pion and 
three-kaon interactions were determined from these calculations. 

5. Conclusions 

The use of lattice field theory methods to compute properties of hadron interactions has rapidly 
progressed in the last few years [25] . Along the path to computing two-baryon interactions, we have 
seen a precision calculation of the 1 = 2 Tin scattering length as well as the computation of systems 
containing up to 12 pions or kaons, allowing the first prediction for both the three-pion and three- 
kaon interactions. In the next few years we will continue to see this rapid progression as the lattice 
community moves into the era of peta-f lop computations. 
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